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We present the high energy neutrino Monte Carlo event generator ANIS 
(All Neutrino Interaction Simulation). The aim of the program is to provide a 
detailed and flexible neutrino event simulation for high energy neutrino detectors, 
such as AMANDA and ICECUBE. It generates neutrinos of any flavor according 
O ■ to a speciflc flux, propagates them through the Earth and in a flnal step simulates 



neutrino interactions within a specifled volume. All relevant standard model 
processes are implemented. We discuss strength and limitations of the program, 
and provide as an example event rates for atmospheric and neutrino spectra. 

1. Introduction 



S ■ The era of dedicated high-energy neutrino telescopes has just begun. The 

^ ■ currently operating detectors AMANDA and BAIKAL will hopefully soon be 

■ joined by ANTARES, ICECUBE, NEMO and NESTOR. The unifying goal is to 



detect neutrinos of extraterrestrial origin. For the Monte Carlo (MC) simulation 
of atmospheric muons in such detectors, CORSIKA is frequently used. The sub- 
\ sequent muon propagation may be done using programs such as MMC, MUM or 

MUSIC. However, no standard for neutrino generators has yet emerged. 

In this paper, the MC z/-event generator ANIS (All Neutrino Interaction 
Simulation) is described. It generates u-events of all flavors, propagates them 
through the earth and flnally simulates i/-interactions within a specifled volume 
around the detector. All relevant SM processes, i.e. CC and NC i/A-interactions 
as well as resonant z/ee~-scattering are implemented. Neutrino regeneration as 
expected in NC-scattering, u + N u + X, and in r production and decay 
chains, z/^ + A ^ r + A, r ^ z/,- + (ui) + X, are included in ANIS. 

In the next sections a description of the characteristic features of ANIS 
is presented. As an application of ANIS, the event rates for atmospheric and 
astrophysical electron neutrinos (following an spectrum) are calculated. 
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2. The Signal Generator ANIS 

ANIS is written in C++ and makes use of the CLHEP package [1]. The 
internal event record is derived from the CLHEP HepMC event class. There 
are currently two output formats implemented: i) the AMANDA specific event 
format f2000 [2] and ii) the HepMC ascii format [1] . 

The program is controlled through a steering file which allows to enable 
specific interaction processes included in the simulation. The currently imple- 
mented interaction channels include CC, NC as well as resonant i^ee~-scattering. 
The cross-section data for CC and NC reactions are provided through pre-calculated 
external tables. The total cross-section is obtained through interpolation. Large 
sets of possible final states, characterized by the Bjorken variable x and y = 
Eh/Ey, have been generated and are randomly sampled from during generation 
of neutrino events. The use of pre-calculated tables makes the program fast and 
independent of other packages. The cross-section data has been calculated up to 
10^^ GeV and is discussed in more detail in the following section. 

Primary i^'s are randomly generated on the surface of the earth with energy 
spectrum F^{E) cx E~^ and are then propagated to the detector. Interacting with 
matter they are either absorbed (CC case) or regenerated at lower energies (NC 
case) [3]. In the special case of CC i^^A^-interaction a short-living r-lepton is 
produced. It propagates in matter, loosing part of its energy, and finally decays 
giving rise to secondary Uj. and, in ~ 17 % of the cases, to secondary or 
z/g. The r-decay in ANIS has been simulated using TAUOLA [4]. Thereby the 
practically full polarization of high energy r's was taken into account. Again, the 
previously generated decays of polarized r-leptons are stored in tables, which are 
then sampled from by ANIS. 

The density profile of the earth in ANIS is chosen according to the Pre- 
liminary Earth Model [5]. Regeneration effects are naturally accounted for by 
feeding the secondary z/'s back into the event record. These z/'s are assumed to be 
emitted parallel to the direction of the primary u. This is justified, since regener- 
ation effects become significant just at very high energies, where the accumulated 
deflection angle is generally smaller than the telescope angular resolution. 

Once the detection volume is reached, a final vertex is sampled along the 
u trajectory within the detection volume (specified through the steering file). In 
the case of a CC t'^A^-interaction, ANIS correctly simulates the muon scattering 
angle. Along with the full event, three weights for later use are written out: the 
normalization of the flux, a weight for the atmospheric flux [6], as well as a weight 
proportional to the total cross-section of the neutrino interaction. 

Event rates for atmospheric and various extraterrestrial neutrino spectra 
are obtained by applying the appropriate weights to the events. This last step is 
to be done by a user defined energy dependent weight function during analysis of 
the events in PAW, ROOT or any other analysis programs. 
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3. Neutrino Cross Sections 

At El, < 10 PeV, deep inelastic z/7V-cross-sections may be successfully de- 
scribed in the framework of pQCD. Parameterization of the uN structure func- 
tions, F^^ [x, Q^), may be taken e.g. according to CTEQ5 [7]. At higher energies 
these cross-sections are dominated by scattering off small x quarks. The unknown 
behavior of the structure functions at x < 10"^ makes the calculations model de- 
pendent. The uncertainty in extrapolations of F^^ [x^Q'^) to small x and large 
influence the expected cross-sections at high energies. 

There are two possibilities presently included in ANIS: i) the smooth both 
in X and in power-law extrapolation of the pQCD CTEQ5 parameterization 
to small X and large Q^, and ii) hard pomeron [8] enhanced extrapolation [9]. 
The cross-sections of the first case, denoted in the Fig. la as pQCD, practically 
coincide with [10], while the second model (HP, dash-dotted curves) predicts 
approximately 2 times higher cross-sections at ~ 1 ZeV. The ^/-distributions, 
normalized to the corresponding cross-sections, are plotted in Fig. lb; the solid 
and dashed curves stand for CC and NC z/iV-scattering, respectively. 

The user can choose the desired model through the steering file. It is 
possible to add any new model to ANIS by providing the corresponding cross- 
section and final state tables. 

Due to rUe <^ ruN, at high-energies practically all i^e~-cross-sections are 
negligible, with the only exception of i/ge" — *■ W~ anything at Ep^ ^ 6.3 PeV 
[11] (see Fig. la). This resonance process is also included in ANIS. Finally it 
should be noted, that ANIS can easily be extended to include new processes 
beyond the Standard Model. 
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4. Discussion 

As an application of ANIS, the resulting distributions of 10^ z/g + v^. events 
for an atmospheric spectrum and a diffuse AGN-like flux [FyiE] = 1 x 10~^£'~^ 
GeV~^ cm~^ s~^ sr~^) are shown in Fig. 2. In order to illustrate the effects of v 
propagation through the earth, the contributions from the upper and lower hemi- 
sphere are shown separately. Alternatively, and can be simulated, however, 
muons produced in the reaction have to be propagated with a separate program, 
e.g. with those mentioned in the introduction. Concluding one can say, that the 
present version of ANIS allows the precise simulation of high-energy i/-events of 
all flavors in a wide energy range and designed flexible to allow incorporation of 
possible new physics processes. 
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